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ABSTRACT 

Aims. The cluster baryon fraction is estimated from the CMB-scattering leptonic component of the intracluster medium (ICM); 
however, the observed cluster baryon fraction is less than the cosmic one. Understanding the origin of this discrepancy is necessary 
for correctly describing the structure of the ICM. 

Methods. We estimate the baryonic mass in the outskirts of galaxy clusters which is difficult to observe because of low electron 
temperature and density in these regions. 

Results. The time scale for the electrons and protons to reach equipartition in the outskirts is longer than the cluster age. Since thermal 
equilibrium is not achieved, a significant fraction of the ICM baryons may be hidden in shells around galaxy clusters. We derive the 
necessary condition on the cluster mass for the concealment of missing baryons in an outer baryon shell and show that this condition 
is fulfilled because cluster masses are comparable to the estimated characteristic mass M = e* Km^^G^) = 1.3x10'^^ solar masses. The 
existence of extreme-ultraviolet emission haloes around galaxy clusters is predicted. 
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1. Introduction 

Clusters of galaxies are the largest virialized structures in the 
Universe. They typically contain hundreds of galaxies, spread 
over a region whose size is roughly 3 Mpc. Their total masses 
are about IO'^Mq. The space between galaxies in clusters is 

filled with very hot ~ 10^ K, low-density ~ 10"^ cm"^ gas (for a 
review, e.g. Sarazin 1986). 

It is widely believed that clusters should be representative 
of baryonic and nonbaryonic matter compositions, and thus that 
their baryon fractions can be used to determine the average den- 
sity of the Universe in conjunction with Big Bang nucleosynthe- 
sis (e.g. White et al. 1993). 

Recent good quality observations of massive clusters with 
Chandra and XMM-Newton put strong observational constraints 
on the gas mass fraction fi, in massive clusters at large radius 
(e.g. Allen et al. 2004, Zhang et al. 2006). In many cases the 
new data allow this fraction to be constrained out to the radius 
rsoo at which the density contrast is 500. The observable baryon 
fraction is lower than the best fit to the WMAP 5-year result of 
ft = 0.17 ± 0.01 (Dunkley et al. 2008). 

X-ray observations are not the only way to constrain the gas 
mass fractions of clusters. The Sunyaev - Zeldovich (SZ) effect 
is also being used for this purpose. LaRoque et al. (2006) esti- 
mated /i, in the radius at which the density contrast is 2500, by 
fitting both BIMA/OVRO SZ effect data and Chandra X-ray data 
simultaneously and fitting the SZ effect data alone. These results 
are in excellent agreement with the X-ray derived results. 

Measurements of the SZ effect in the CMB maps of WMAP 
data (Afshordi et al. 2007), through study of a sample of 193 
massive galaxy clusters with observed X-ray temperature greater 
than 3 keV, indicate that a significant fraction, 35% + 8%, of 
baryonic mass is missing from the hot ICM, and thus must cool 
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to form galaxies, intracluster stars, or an unknown cold phase of 
the ICM. There does not seem to be enough mass in the form 
of stars or cold gas in the cluster galaxies or intracluster space, 
signaling the need for a still unknown baryonic component, or 
otherwise new astrophysical processes in the ICM (Afshordi et 
al. 2007). 

Much larger discrepancies with X-ray observations have 
been reported in Lieu et al. (2006), who claim that the WMAP 
SZ signal for ROSAT X-ray clusters is a factor of 4 weaker than 
the expectation from X-ray brightness profiles. 

However, in cosmological simulations the baryon fraction 
within rsoo has converged to 90-95% of the best fit to the WMAP 
5-year data (Dunkley et al. 2008) in the case of non-radiative 
simulations (e.g. Frenk et al. 1999, Kay et al. 2004, Crain et al. 
2006). Cosmological cluster simulations that include gas cooling 
and star formation also give an ICM mass fraction very similar 
to the WMAP 5-year fraction (Kravtsov et al. 2005). 

In a hierarchical Universe clusters form through accretion. 
As fresh material joins the stationary cluster gas, it converts its 
infall kinetic energy into heat through a shock. A shock front 
moves outward as ambient gas accretes toward the cluster cen- 
ter (see e.g. Takizawa & Mineshige 1998). The accretion shock 
around an X-ray cluster primarily heats the protons, since they 
carry most of the kinetic energy of the infalling gas. The elec- 
trons, on the other hand, remain cold after the shock with their 
temperature slowly rising as a result of their Coulomb collisions 
with the hot ions. Thermal equilibrium will only be achieved if 
the temperature equilibrium time between the protons and the 
electrons is shorter than the age of the cluster (e.g. Fox & Loeb 
1997). 

The Sunyaev-Zeldovich distortion of the CMB spectrum is 
proportional to the electron temperature. Thus, the low gas frac- 
tions inferred from the SZ effect observations in the CMB maps 
of WMAP data (Lieu et al. 2006, Afshordi et al. 2007) can be 
interpreted as evidence of a warm phase (with low electron tem- 



2 



D. A. Prokhorov: Missing baryons in shells around galaxy clusters 



perature Tg ~ 0.1 keV and high proton temperature) of the ICM 
in the cluster outskirts, which would not contribute appreciably 
to the SZ signal, while containing a significant fraction of the 
ICM baryons. 

The importance of cluster outskirt analysis is discussed in 
Sect. 2, where the baryon fraction profile is calculated in the 
case of a simple isothermal beta model. The formation of bary- 
onic shells around clusters is considered within the framework 
of kinetic theory in Sect. 3. A condition on the cluster mass for 
deviation from thermal equilibrium is derived in Sect. 4. The 
influence of the deviation from thermal equilibrium on the pre- 
diction for the SZ decrement is considered in Sect. 5. A method 
for detecting missing baryons by means of EUV astronomy is 
proposed in Sect. 6 and a discussion is presented in Sect. 7. 

2. Tracing the baryon fraction in the outskirts within 
the framework of the hydrostatic theory 

The ICM emits energy mainly via thermal bremsstrahlung, 
which is proportional to the density squared. Therefore external 
cluster regions show less X-ray emission than the central ones, 
where the density is much higher This lower emission translates 
in lower statistics for available X-ray observations. One possible 
way to overcome the problem of low statistics in the outer re- 
gions is to fit clusters with the isothermal beta model (Cavaliere 
& Fusco-Femiano 1976). The disadvantage of such a method 
is clear: the extrapolation depends on the validity of the model 
in the cluster outskirts. In accordance with the beta model, the 
electron number density follows the spherical distribution 



n{r) = «o 1 + — 



(1) 



where no is the central number density, Vc is the core radius, yS is 
the slope parameter 

The hydrostatic equilibrium equation for a spherically sym- 
metric isothermal cluster (e.g. Sarazin 1986) can be written as: 



kTdn — — 



GM(R)iMipn(R)dR 



(2) 



where M{R) — Anr ip{r) dr is the total mass within a sphere 
of radius R, i^(r) is the mass density of all the matter, yU the mean 
molecular weight, nip the proton mass. 

The baryonic fraction /i(r) = p{r)l^{r) is given by 



Mr) 



kT (r^p-Hr)- dp(r)/dr)' 



(3) 



where p(r) - mpn(r) is the gas mass density, the prime denotes 
a derivation of r 

Using Eq.([TJ, one can calculate the ratio fb{r)/fb{rc) for the 
isothermal beta model 



fbirc) 



2\2-3/?/2 



X 3 + 



(4) 



The maximum of the function F{r) is at the point a = r/r^: 



2/3 <^< 10/9. 

The maximum value of F{r) is given by 
Fmax^ -y/(4-3y6)4-3/?(3/3- 2)3/^-2. 



(5) 



(6) 




Fig. 1. Position of the maximum of the fbif)/fb{rc) ratio as a 
function of the slope parameter j8. 




Fig. 2. Maximum value of the fb{r)l fbirc) ratio as a function of 
the slope parameter p. 



The dependence a(y6) and the maximum value of Fmax(J3) 
are plotted in Figs. 1 and 2 respectively. 

If yS < 2/3 then the baryon fraction F{r) is an increasing 
function of the radius. The average value of /3 found by Jones & 
Forman (1984) from the X-ray distribution is yS = 0.65, which 
implies that the gas is considerably more extensively distributed 
than the mass in the cluster and the baryon fraction in the out- 
skirts is higher than in the central region. Hence a study of the 
cluster outskirts is essential for understanding the baryon com- 
ponent distribution in clusters. 

The mass of the gas shell between the accretion shock radius 
around a galaxy cluster and the virial radius 7?,,,,. in the context of 
the beta model (with /? < 1 and «: /?,,,>) is given by integrating 
Eq.(l): 



Mf;;:„.M^-(R, 



(7) 



where M^"\Ryir) is the gas mass inside the virial radius sphere, 
and a is the ratio of the shock radius to the virial radius. 

Lieu et al. (2006) claim that the WMAP SZ signal for 
ROS AT X-ray clusters is much weaker than the expectation from 
the X-ray brightness profile. In order to compare with WMAP 
maps, they extrapolate the isothermal beta model fits far beyond 
the region that is fit by X-ray data. The reported discrepancy may 
be evidence of a warm phase of the ICM in the cluster outskirts 
with electron temperature ~ 0.1 keV, which would not con- 
tribute appreciably to the SZ signal, and the proton temperature 
is of the order of the virial one (see Sect. 5). 
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3. Formation of a baryonic shell within the 
framework of the kinetic theory 

In accordance with the hierarchical clustering scenario, the 
baryons are accreted along with the dark matter onto a galaxy 
cluster An accretion shock at the interface of the inner hydro- 
static gas with a cooler external medium is a long-standing pre- 
diction from gravitationally driven models (e.g., see the model 
of Bertschinger 1985). 

The accretion shock converts the radial motion energy of the 
upstream protons into random motion energy, below the shock. 
Let us assume that the protons below the accretion shock are es- 
sentially collisionless and have large orbital eccentricities. Then 
they are able to accrete into the X-ray cluster on the Coulomb 
time scale, since this is the time in which they are able to dissi- 
pate their excess orbital angular momentum (for a similar view 
of a quasar model, see Meszaros & Ostriker 1983). If this time 
scale is comparable with the cluster age, then a considerable part 
of the protons can be located between the accretion shock and 
the inner collisional plasma region. The Coulomb mean-free- 
path of protons is larger than the virial radius of the cluster at 
temperature > 2 keV (Loeb 2007), therefore it is reasonable to 
choose the virial radius as the inner collisional region border. 

The Coulomb mean-free-path of a proton in a plasma of tem- 
perature T is given by (Spitzer 1956): 



An^lhie^ In A 



(8) 



where T is the proton temperature, n is the proton number den- 
sity. In A is the Coulomb logarithm, e is the electron charge, ks 
is the Boltzmann constant. 

The proton random energy (or "temperature") is of the order 
of the gravitational temperature. 



T - 

1 gr — 



GM{R,i,.)m„ 

k-BRvir 



(9) 



where M{Ryii) is the total mass inside the sphere of radius /?,,,,., 
which equals (e.g. Eke et al. 1998) 



Rvir 



\ 4-nAcPcr 



1/3 



(10) 



where A^. is the density contrast for the formation of the cluster. 

The protons will be collisionless if the Coulomb mean-free 
path / is higher than the distance L, which the proton can cover 
during the shock propagation time t - tH^^ from the virial ra- 
dius to the present accretion shock radius: 



/ > L 



(11) 



where Hq is the Hubble constant. 

It is convenient to rewrite the distance L by means of Eqs. 
^ and ((TOl): 



The mass of the gas shell between the accretion shock and 
virial radius is 



M^-;, = —Ri 

she II 2 I' 



(3_ 

R\ 



- 1 



(14) 



where R^ is the present shock radius. 

The overdensity in Bertschinger's collisionless model, given 
by the dimensionless quantity M{X)ly\?, equals 178 at /l„> ^ 
0.255 (Fox & Loeb 1997), where the dimensionless parame- 
ter A - rjr,a and r,a is the turnaround radius. The shock oc- 
curs at fixed /I, ■x, 0.347 (for j - 5/3) in accordance with 
the shocked accretion model of a collisional gas (Bertschinger 
1985). Consequently, the ratio of the shock radius to the virial 
radius a - RJRyi,- equals 1.36. 

One can rewrite the condition / > L in terms of the shell mass 



shell 1, \ ' ^" ' I 



replacing T by Tgr from Eq. (|9]l one can get: 



shell 



< T 




1)' 



, G'-M\R,i,) 



(15) 



(16) 



If the missing baryons are situated in the collisionless shell, one 
can evaluate the value M7,,/M*°^(7;,,,>) = {fi, - f)lf, where 
/ = M«"\R,ir)IM(R,ir) and f,, is the best fit to the WMAP 5- 
year result. Therefore Mf^-^„/MiR,i,) = fb-f. 

Thus a collisionless shell should exist if the total cluster mass 
satisfies the necessary condition 



M{R,i,)>T- 



lA,{ fh-f)lnA / 
- 1 G^n 



(17) 



The time evolution of the radius rs(t) of a shock surface fol- 
lows rs{t) oc fl'^ (Bertschinger 1985). The dimensionless shock 
propagation time from the virial radius to the present accretion 
shock radius equals t = 1 - a^'^l^, assuming here that the cluster 
age is the Hubble time. 

If fb - f - 0.056 (if 33% of baryons are missing), a - 1.36 
(the ratio of the shock radius to the virial radius in Bertschinger's 
model) and r = 0.29 (the shock propagation time), then the total 
cluster mass must exceed the limit mass value 



M{Ryir) > 1.2- IO'^Mq. 



(18) 



X-ray measurements tend to come from the inner parts of 
clusters. One way to relate an observable quantity in the inner 
regions of a cluster to a quantity in the outer regions is to use a 
simplified model, such as that considered by Eke et al (1996), in 
which both the X-ray gas and the total mass have the form of a 
singular isothermal sphere in the virialized section of the cluster 
In this model, the mass within the sphere of the virial radius and 
the temperature of the X-ray gas are related by 




(12) M{R,,,)^ 1.7 



10'% 



3/2 



M, 







(19) 



Using Eqs. (O and (fTTT i. one can show that the collisionless 
condition is equivalent to the inequality 



n < 



(13) 



where h-joi is the Hubble constant in units of 70 km- s ' Mpc ' . 

The limit mass value (1.2 ■ IO'^Mq) is comparable with 
a typical cluster mass, therefore a fraction of baryons can be 
trapped in the baryonic shell around clusters within the frame- 
work of kinetic theory. 
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4. Deviation from thermal equilibrium 

In Sections 2 and 3 using both hydrostatic and kinetic theories 
we conclude that the fraction of cluster baryons located in the 
outskirts may be significant. This raises the possibility that a 
missing fraction of baryonic mass may be hidden in the cluster 
outskirts where the electron temperature is much smaller than 
the proton one, and therefore the contribution of these electrons 
to the SZ signal is small. 

The time scale teg{p, e) for the electrons and protons to reach 
equipartition - Tp is given by teqip, e) k yjniplme ■ t^qip, p) 

(e.g. Sarazin 1986), where the time scale teqip, p) - U s/kTgr/nip 
for protons to equilibrate among themselves . 

The electrons will be cold between the accretion shock and 
virial radius if the time scale t^qip, e) is longer than the shock 
propagation time t - t//q ' from the virial radius to the present 
accretion shock radius, i.e. t^qip, e) > tH^^ . This condition can 



be rewritten in the form equivalent to Eq.lflll) 




(20) 



Assuming that the missing baryons are situated in the shell 
between the accretion shock and virial radii and following the 
same method as in Section 3, one can derive a necessary condi- 
tion for deviation from thermal equilibrium: 



/)lnA e* 




(21) 



fh - f - 0.056 (if 33% of baryons are missing), a - 1.36 
(the ratio of the shock radius to the virial radius in Bertschinger's 
model) and t - 0.29 (the shock propagation time), then the nu- 
merical value of this mass is 



M{R,ir) > 2.8 ■ lO'^M, 



o- 



(22) 



The value of the limit mass M(R,ir) = 2.8-10"Mq is smaller 
than cluster masses for a sample that was used by Afshordi et 
al. (2007) and therefore the missing baryons may be hidden in 
shells around galaxy clusters. Note that even for a value t = 1 
(the cluster age) the limit mass is still smaller than cluster masses 
for the Afshordi et al. (2007) sample. 

This explanation of the origin of the missing baryonic mass 
in the clusters arises from a coincidence between a typical cluster 
mass and a characteristic mass (1.3 ■ 10''' Mq), which is given 
by 



M = 



mlG^' 



(23) 



The characteristic mass depends only on fundamental phys- 
ical constants. 



5. Prediction for the SZ decrement 

The observed SZ effect with WMAP is much less than the pre- 
dicted decrement value given by the isothermal beta model (Lieu 
et al. 2006). Their study used a single-temperature model for the 
intracluster medium, under the assumption that equipartition is 
fully achieved. However, the outskirts of galaxy clusters actu- 
ally have a two-temperature structure, because the two tempera- 
tures, i.e. those of electrons and protons Tp, remain different 



(e.g. Fox & Loeb 1997). To understand the discrepancy between 
observed and predicted SZ decrements it is useful to study the 
influence of the deviation from thermal equilibrium on the pre- 
diction for the SZ decrement. 

The thermal SZ effect (see Sunyaev & Zeldovich 1980 for a 
review) describes the inverse Compton scattering of CMB pho- 
tons by the hot gas in clusters. The resulting change in the (ther- 
modynamic) CMB temperature observed at frequency v is given 
by 



l^Ty = yj(x)To 



(24) 



where Tq = 2.725K is the unperturbed CMB temperature, y is 
the Comptonization parameter, x is a dimensionless parameter 
defined as x = hv/kTo and in the nonrelativistic regime {kTe « 



+ 1 

- 1 



■4. 



The Comptonization parameter is 

dl — -c 

J mc~ 



y 



ro-jn 



(25) 



(26) 



where is the electron temperature, ctt is the Thomson cross 
section, and the integral is over the distance along the line-of- 
sight. 

The total flux Sy of a. cluster (e.g. Refregier et al. 2000) ob- 
served at frequency v is related to the angular integral of the 
thermodynamic temperature shift hTy by 



Sy = -^^q{x) I dD.My 



(27) 



where qix) = x4/[2 sinh(x/2)]2. For the WMAP W channel 
(v - 94GHz, X94 - 1.65): ^(^94) = 2.18 and the fufl width at 
half maximum (FWHM) equals 10 arcmin. 

The postshock protons carry nearly all of the bulk kinetic 
energy. After passing through a shock the electrons can gain 
thermal energy via collisions with protons. The evolution of the 
electron temperature is given by 



dTe _ Tp-Te 
dt teqip, e)' 

The collision time scale t^qip, e) (Spitzer 1958) equals 



teqip, e) = 



3memp 



8 V27r«e4jjjyv\ m, 



+ 



3/2 



(28) 



(29) 



Using the dependence of the shock radius r^it) cc r^^^ from 
the work of Bertschinger (1985) and Eqs. ( |28] |. (l29l l and assum- 
ing typical values of the central electron temperature ksTo - 8 
keV, of the central density «o = 10"^ cm"^, and of the slope 
parameter /? = 2/3, we found the electron temperature profile 
and calculated both the SZ temperature decrement and the total 
SZ flux Sy within 10' radius for the non-isothermal beta model. 
For a cluster with a core radius 0.5' which corresponds to the 
mean core radius in the sample of clusters by Bonamente et al. 
(2006), the SZ decrements are plotted in Fig. 3 for the isothermal 
and non-isothermal beta models. The total SZ flux of a cluster in 
the non-isothermal model is smaller by a factor of 2.6 than those 
computed with a single-temperature model. The SZ decrement 



' Note that this is the connect formula, while there is a misprint in the 
expression of q(x) in Eq.(19) by Refregier et al. (2000). 
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X [arcmin] 

Fig. 3. SZ temperature decrements: for the isothermal beta 
model (solid line), for the non-isothermal beta model (dashed 
line), and for the isothermal beta model in which the value of 
the total SZ flux corresponds to the non-isothermal model (dot- 
dashed line) 



for the isothermal beta model which is scaled in such a way that 
the total SZ flux corresponds to the non-isothermal model is also 
plotted in Fig. 3. 

Although the total SZ flux for the non-isothermal beta model 
is significantly lower than that for the isothermal model, atten- 
tion should be drawn to a paper of Bielby & Shanks (2007), 
where they truncated the isothermal beta model at the 2' radius 
and still found a predicted SZ effect much larger than the level 
observed by WMAP, for a sample of clusters (Bonamente et al. 
2006). Note that the approach of Afshordi et al. (2006) is dif- 
ferent from that of Bielby & Shanks (2007), in the former the 
X-ray data are mainly used to define a template to detect SZ 
decrements and then the SZ data and the X-ray temperature data 
alone are used to establish the gas densities. 



6. EUV emission from a baryonic shell 

Extreme ultraviolet (EUV) astronomy has produced substantial 
results in a wide variety of fields (e.g. Bowyer et al. 2000). 
The Extreme Ultraviolet Explorer (EUVE) has detected emis- 
sion from a few clusters of galaxies in the 70-200 eV energy 
range. 

Since baryonic shells around clusters are a possible location 
of missing cluster baryons and since the time scale of collisions 
between baryons and electrons is longer than the Hubble time in 
such shells, these cold electrons should be a source of extreme- 
ultraviolet emission. 

In this section the contribution to the total cluster EUV emis- 
sion from a baryonic shell is estimated for the Coma cluster The 
Coma cluster is a rich, hot (r,„=8 keV), nearby {D - lOO/i^^ 
Mpc) galaxy cluster, which is often considered as a "standard 
cluster" due to the wealth of observational data available for this 
object. 

The spectral surface brightness along a particular line of 
sight writes as 



b{E) 



(r)A(£, T,)dl 



(30) 



where A(£', Tg) is the spectral emissivity of the gas at observed 
energy E, and is the electron temperature. 



The spectral emissivity due to the electron-proton 
bremsstrahlung rate is given by 



ME, T,) 



1 



exp 



(31) 



For the sake of clarity, we propose that the hot gas is lo- 
cated inside the virial sphere. The electrons remain cold inside 
the baryonic shell, because the time scale of their collisions with 
the hot protons is longer than the Hubble time. The temperature 
of these electrons is almost the same as the preshock tempera- 
ture, a plausible value is Text = 3 ■ 10^ K. 

The gas density is given by the isothermal beta model (Eq. 1) 
within the virial sphere, and the gas density inside the shell, 
which corresponds to missing baryons in the galaxy cluster (if 
33% of baryons are missing), equals 



nshell 



: 0.5 f « 

Jo 



47T 



4-7T 



ir)4nr^dr-\-Ri--R:.^ 



(32) 



where the center number density no = 3.4 ■ 10"^ cm"^, the core 
radius r^. = 290 kpc, /3 - 0.75, the virial radius /?,,,, = 2.8 Mpc 
for a Coma-like cluster (Hq - 70 km- s"' Mpc"' is adopted). 
The spectral surface brightness from the hot gas is 



bi,{0) 



1 + 



sin- OV''' e-^'^"- D 



r? cos^ O 



Tin COS- <S> 



(33) 



where is the angle between the center of the cluster and the 
direction of interest <1>„m.i is defined by 



arctg 



DsinO 



(34) 



where D is the distance of the cluster 

The spectral surface brightness from the cold electrons is 
given by 

if < Ryir/D: 



bc{0) ~ 2 



•J<s>„„,Ae) 



cos^ <1> 



if 6* > Ryi, /D: 



(35) 



(36) 



where (S>ma.a is 



'^maxi = arctg 



^Rl - D- sin^ 



DsinO 



(37) 



B(0) = 



The normalized spectral surface brightness writes as 
bM + bM 



b,M + bM' 



(38) 



For the Coma cluster, the normalized spectral brightness at 
observed energy 0. 1 keV is plotted in Fig. 3, where 1° is 1 .7 Mpc. 
Fig. 3 shows the presence of an extreme-ultraviolet emission halo 
around the galaxy cluster without the X-ray emission from it, 
because only the cold electrons are situated in the baryonic shell. 
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Angle (degree) 

Fig. 4. Normalized surface brightness for the Coma cluster: the 
hot gas (solid line), the hot gas + the baryonic shell (dashed line). 



The cold electrons in shells around galaxy clusters should be 
a source of extreme-ultraviolet emission. Observations of cluster 
outskirts can shed light on a location of missing baryons. 

Other possibilities that have been proposed to solve the clus- 
ter missing baryon problem are; diffuse intracluster light, which 
would be challenging to detect observationally, but could poten- 
tially hide significant amounts of stellar mass in the vast intra- 
cluster space (Lin & Mohr 2004, Gonzales et al. 2005, Zibetti 
et al. 2005, Monaco et al 2006); hiding baryons in the form of 
cold and compact dark baryonic clouds from local cooling insta- 
bilities within the ICM (Dwarakanath et al. 1994); evaporation 
of baryons out of the virial radius of the cluster (Loeb, 2007); 
the correct value of Q.m is higher than the best-fit WMAP value 
(McCarthy et al. 2007). 

Acknowledgements. I am grateful to Florence Durret and Roya Colin Mohayaee 
for valuable discussions and to the referee for constructive comments. 



7. Discussion and conclusions 

Cluster gas mass fractions as inferred from X-ray observations 
have been used as a probe of the universal ratio of baryon to total 
mass densities Q.bI^m (e.g. White et al. 1993, David et al. 1995, 
Evrard 1997, Roussel et al 2000, Allen et al. 2002, Ettori 2003, 
Allen et al. 2004). But the result of recent measurements of SZ 
effect and X-ray studies is that a significant fraction 35% ±8% of 
the baryons is missing from the ICM (Ettori 2003, LaRoque et 
al. 2006, Vikhlinin et al. 2006, Afshordi et al. 2007, McCarthy et 
al. 2007). A possible explanation is that a significant fraction of 
baryonic mass may be hidden in a form that is difficult to detect 
with ordinary observational means. 

The problem is unlikely to be associated with the conven- 
tional cluster mass estimates at least within the Abell radius, 
which rely upon the hydrostatic equilibrium hypothesis for the 
dynamical state of clusters in the computation of their total 
masses (e.g. Wu & Xue 1999). This point has been justified by 
the excellent agreement among the X-ray, optical and weak lens- 
ing determined cluster masses on Abell radius scales (e.g. Allen 
1998, Wu et al. 1998). However how accurate the conventional 
cluster mass estimates will be in the outskirts of clusters where 
hydrostatic equilibrium becomes questionable remains unclear. 

We conclude within the framework of hydrostatic theory that 
the gas is considerably more extensively distributed than the 
mass in the cluster. Therefore a significant fraction of the ICM 
baryons may be located in shells around galaxy clusters. 

In this paper the formation of baryonic shells around clus- 
ters is considered within the framework of the kinetic theory. 
Baryonic shells can be formed in this way for galaxy clusters 
with masses (M > 1.2 ■ IO'^Mq). 

We propose that the low baryon fraction inferred from the SZ 
effect observations in the CMB maps of WMAP data (Afshordi 
et al. 2007) can be explained by the deviation from thermal equi- 
librium in the cluster outskirts. The electrons and protons do not 
reach equipartition there. This warm phase (with electron tem- 
perature ~ 0.1 keV and with high proton temperature) of the 
ICM in the cluster outskirts does not contribute appreciably to 
the SZ signal. The necessary condition on the cluster mass for 
the concealment of the cluster missing baryons in a baryon shell 
is fulfilled because cluster masses are comparable to the value of 
the characteristic mass M - e'^/{m^pG^) = 1.3 • IO'^Mq. 

We considered a non-isothermal beta model for the ionized 
gas distribution of a cluster and showed that the total SZ flux 
of a cluster differs by a factor of 2.6 from that computed with a 
single-temperature beta model. 
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